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President’s Corner 


By Steven Biste, N7HPR, Vice Presipent, TAPR 


T’m writing the President’s 
Corner this quarter while Dave 
is recuperating from surgery. 
For those of you who know 


Dave, you know he won’t be down long. We all wish 
him a speedy recovery. He did say he would be at 
Hamvention. See, nothing can keep him down! 


invol 


kit is 


availal 


at have we been up to this past quarter? TAPR’s 
lvement with HPSDR remains strong. The LPU 


ble and can be ordered online. Sourcing 


parts for Alex has been challenging. We have been 
working with an off shore contract manufacturer to 


wind 
sure 
time 
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ly read 
ay at tl 


Loo 


the many coils. It’s 
they are done correctly. This has taken more 
than we thought. The Pandora enclosure is 


a tedious process to make 


ly for production and prototypes will be on 
¢ TAPR Hamvention booths (#455458). 


for upcoming production announcements for 


the Pennywhistle 20W PA and the Excalibur 10MHz 
Frequency Reference. 

The new TADD-2 frequency divider kit is now 
available. It accepts a 5 or 10 MHz input and provides 
six low impedance outputs. Each can be set to a rate 
of 1 PPS (with two duty cycles), 10 PPS, 100 PPS, 
1000 PPS, or 10K PPS. For more information, go to 
page 3. 

‘What are we up to next quarter? Hamvention is 
only days away. We have a great line up for the Friday 
morning TAPR Digital Forum starting at 9:15 AM in 
Room 1. 

* Intro and TAPR Update - Steve Bible, N7HPR 


* Timing Products Update - John Ackermann, 
N8UR 

¢ HPSDR Manufacturing Update - Scotty 
Cowling, WA2DFI 


* From Napkin to Ham Shack: Creating a TAPR 
Project - John Ackermann, N8SUR 
* Sharing Your Experience, Documenting Your 
Journey - Larry Wolfgang, WR1IB 
* ARRL/TAPR Digital Communications 
Conference Update - Steve Bible, N7HPR 
‘We are pleased to announce that Richard Garriott, 
W5KWAQ, will be the keynote speaker at the 
AMSAT/TAPR Banquet Friday night. Be sure and 
get your tickets in advance. The deadline is May 11th. 
Tickets will not be available at Hamvention. To order, 
see the announcement later in this newsletter. 
Come see us at the TAPR booth and see what 
we are up to. Join in the fun, talk and enjoy the 
camaraderie. 
~ Steve, NTHPR 
eee 
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TAPR at Hamvention 


By Stan Horzerpa, WA1LOU 


TAPR will be at Hamvention in Dayton, Ohio, 
May 15-17. 
TAPR’s Bootn Space 

See TAPR at booths 455-458 (next to AMSAT’s 
booths) in the Ball Arena, where we will have 
HPSDR equipment up and running along with 
other TAPR goodies, like the new TADD-2 
frequency divider kit. 
TAPR Dicitat Forum 

The TAPR Digital Forum will be 9:15 to 11:15 AM 
on Friday, May 15 in Room 1. Steve Bible, N7HPR, 


will moderate and the presentations will include: 


ntroduction and TAPR Update” by Steve Bible, 
N7HP! 


“Timing Products Update” by John Ackermann, 
N8UR 


“HPSDR Manufacturing Update” by Scott 
Cowling, WA2DFI 


“From Napkin to Ham Shack: Creating a TAPR 
Project” by John Ackermann, NSUR 


“Sharing Your Experience, Documenting Your 


Journey” by Larry Wolfgang, WR1IB 
“ARRL/TAPR Digital Communications 


Conference Update” by Steve Bible, NTHPR 


Immediately following, TAPR Secretary, Stan 
Horzepa, WAILOU, moderates the APRS Forum 
at 11:30 AM to 12:30 PM, also in Room 1. 


TAPR/AMSAT Banauet 


he third annual TAPR/AMSAT Banquet will be 
on Friday evening, May 15, 2008. 


he banquet will be in the Jefferson Room at 
the Kohler Presidential Banquet Center, 4572 
Presidential Way, Kettering, OH (just south of 
Dayton). 


Richard Garriott, W5KWAQ, will be the keynote 
speaker. Richard was aboard the ISS (International 
Space Station) in October 2008 and delighted many 
Amateur Radio operators and school children with 
contacts. He was the sixth private citizen to make 
the trip. In 1983, Richard’s father Owen Garriott, 
W5LEL, was the first person to make Amateur 
Radio contacts from the Space Shuttle. 

Reservations are required and must be made by 
Monday, May 11, 2009; tickets will not be sold at 
Hamvention. The price for the banquet is $30 per 
person. Reserved tickets may be picked up at the 
AMSAT booth at Hamvention on Friday, or at the 
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door Friday evening. 


Make reservations online in the AMSAT Store 
at http://www.amsat.org or phone Martha at the 
AMSAT office, 10 AM to 6 PM EDT at 301-589- 
6062 or toll free (in the US) at 888-322-6728. 


Banauet DetaiLs 


6:30 PM: Doors open and cash bar is available 
with beer, wine, liquor and soft drinks. 


7:15 PM: Buffet dinner service begins. 
The menu includes: 
* Garden salad with assorted dressings 


* Barbeque Baby Back Ribs Marinated 


* Roasted Garlic Rosemary Chicken Breast in 
lemon butter sauce 


* Salmon with Newburg sauce 


* Petite Double Stuffed Baked Potatoes 


* Normandy blended green beans 


* Fresh fruit bowl 


* Roll and butter 
Coffee/ iced tea/ water 
* Assorted pies 


ane 
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New Time and Frequency Products: 
TADD-2 and TADD-2 Mini Frequency Dividers 


By Jonn AckeRmaNN NSUR 


Frequency counters are great tools. But they have 
some limitations: they aren’t very good at describing 
changes in frequency over time, and in fact they aren’t 
of any use at all in telling what time it is. 


Asa result, many time and frequency measurements 
rely on phase rather than frequency. Phase is just the 
difference in time between two “clocks.” By measuring 
the rate at which one clock ticks versus the other, we 
can learn a lot about their relative performance, and 
that data can be converted back to frequency error if 
desired. 


While phase can be measured at any frequency, it’s 
often very convenient to use one pulse per second 
- abbreviated as a “PPS” signal. For example, GPS 
receivers used for timing generate a very precise PPS 
signal. While some clocks and frequency standards 
provide a PPS output, many don't. It’s not too hard 
to wire together a bunch of divider chips to turn a5 
or 10 MHz signal into a pulse-persecond output, but 
some of the chips commonly used drift quite a bit with 
temperature, one-off IC projects can be clunky, and 


they usually don’t have the right configuration for the 
next experiment you want to do. 


‘When I designed the TADD-1 distribution amplifier, 


TADD-2 Frequency DivinerR 


my intention was to do a PPS divider project next, and 
to give it the flexibility a useful lab instrument should. 
However, I got bogged down in an overly complicated 
design using technology that was new to me (a CPLD 
chip), and at the time I needed a PPS distribution box 
for my existing signal sources more than a divider. So, 
the planned TADD 2 was put on the back burner, and 
the TADD3 came out instead. 

Anumber of years ago, one of the “time-nuts” gang, 
Tom Van Baak, wrote a very clever program for a PIC 
processor that acted as a singlechip frequency divider. 


Clocked with a 10 MHz input, it would generate 
outputs at 1 PPS as well as several other frequencies, 
with very little jitter. It also had the ability to synchronize 
itself to an external pulse. However, Tom’s code was 
limited to a 10 MHz input frequency; while that is 
very useful, a lot of frequency standards have a 55 MHz 
output. Some time later, another timenut, Richard 
McCorkle, modified Tom’s code to work at either 5 
or 10 MHz. Tom and Richard both graciously agreed 
to Open Source their code, and a few months ago I 
started work on the long-delayed TADD22, designed 
around a PIC divider. A follow on from that project, 
the TADD22 Mini (T2-Mini) is also in the works. 


TADD-2* 


TAPR is now pleased to announce that the TADD-2 
frequency divider kit is now available. Member price 
will be $62; non-member price will be $69. 


* In this case, “TADD” stands for “Totally Awesome 
Digital Divider.” As you may recall, for the TADD-1 it stands 
for “Totally Awesome Distribution Device” and for the TADD- 


3 it stands for “Totally Awesome Digital Distribution.” Here at 
TAPR, we are doing our part to preserve our natural acronym 


resources. 
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The TADD22 accepts a 5 or 10 MHz input (jumper 
selectable) and provides six low impedance outputs 
of about 2.75 volts into a 50ohm load. Each output 
can be set to a rate of 1 PPS (with two duty cycles), 10 
PPS, 100 PPS, 1000 PPS, or 10K PPS. The outputs are 
normally low, going high on the leading edge of the 
pulse. However, each can be inverted (i.e., normally 


high, going low on the leading pulse edge). There is also 


aseventh auxiliary PPS output available that is TTL 
level, 1 PPS only, and cannot be inverted. 


The input circuit allows signal levels as low as -10 dBm. 


to clocl 


the divider. Like the TADD-1, the RF input is 
transformer coupled and DC isolated. A jumper allows 
either 


input. 


igh impedance or a 50-ohm termination for the 


A synchronization circuit allows the divider to 
be stopped, reset, and then synchronized with the 
leading edge of an incoming pulse. This allows the 
output pulse to be aligned with, for example, a GPS 
signal. Synchronization is a twostep process. Pressing 
the ARM switch causes the counter to stop, reset, 
and wait for the arrival of the SYNC pulse. After 
synchronization, the output will be within one clock 
oycle (e.g., 100 nanoseconds at 10 MHz) of the SYNC 
pulse. The SYNC pulse must be at TTL levels, but its 
polarity can be jumper selected. 
The TADD22 is the same size (roughly 6 x 4 inches) as 
the TADD-1 and TADD3, and will fit into the TADD 


enclosure. It runs from a 12-volt supply (with onboard 


voltage regulators) and draws from 50 to about 300 MA = The TADD22 uses through-hole components and 


depending on the output load and polarity. 


Preliminary measurements of the TADD-2’s jitter 
(in other words, how much the PPS output wanders 
around due to circuit noise) shows a standard deviation 
of about 3040 picoseconds, not much more than the 
noise of the best available time interval counters. 

I mentioned above that phase shift over temperature 
could be a problem with dividers made out of 
cascades of ICs. I have done some very preliminary 


TADD-2 Mini Frequency Diviver 


measurements on the TADD2, and it seems to 
perform very well in this regard. Measuring from 
room temperature up to about 73 degrees C, the PPS 
output changed phase by about 3.25 nanoseconds, 
or 60 picoseconds per degree. I plan to rerun these 
measurements shortly with improved temperature 
monitoring to get more reliable results. 


will be available as a kit only. A PIC programmed with 
the divider code is included, and the source code is 
available for those who like to tinker. 


TADD-2 Mini 
In talking with Tom Van Baak about the TADD2 


design, we came up with the idea of a smaller version 
of the divider, a “dongle” that could easily be used for 
temporary experiments, or mounted inside an existing 
frequency standard. I decided to create a small, simple 
version of the divider. The end result is a circuit board 
that is 2.0 inches long by 0.75 inches high. It is built 
using surface mount components (except for the PIC 
chip itself, which is an 8 pin DIP package). 

The TADD22 Mini (or “T2-Minj”) is a much simpler 
device than the fullsize TADD 2. It has only a single 
low impedance output channel, and provides only 

a PPS output (no higher frequencies). The output 
polarity isn’t selectable. Changes in input frequency 
require changing solder jumpers rather than moving 
shorting blocks. On the other hand, it’s tiny and it’s 
cute. 


The T2-Mini is still under development, but is 
progressing rapidly. Because it uses small surface mount 
components that may be difficult for some users to 
handle, TAPR will offer it as an assembled and tested 
unit only. The price has not yet been determined. 

TAPR hopes to start shipping the T2-Mini this 


summer. 


ane 


APR PSR 08 SPRING 200 


Treasurer’s Report for 2008 


Susmittep By Tom Hoimes, N8ZM, TREASURER 


For the fiscal year ended 12/31/2008, TAPR’s 
finances look to be in very good shape, but there 
are some caveats to keep in mind. The biggest one is 
that while we appear to have over $135,000 in assets 
that is mainly because of our sales model, which 
takes in sales dollars before we pay for the parts that 
go into the HPSDR boards. In this case, Mercury 
sales just in the month of December totaled almost 
half of our total kit revenue for the year, but we had 
yet to pay for any of the parts or assembly labor that 
goes into it. The evidence for this is that our total 
December expenses for all categories amounted 


to less than onethird of the expenses for the year. 
When the Mercury bills are all paid, the books will 
show a more realistic picture of our finances. 


This is not to imply that we are in dire straights, 
but merely to point out that the timing of this 
product distorts the picture a bit. A more normal 
picture of our situation would show us with about 


$40,000 combined from our checking and savings 
accounts, and the P&L situation to be enough 
better than break even to allow us to fund future 
development projects comfortably. 


Given that the overall picture is healthy (kind of 
unusual in the current business climate, don’t you 
think?), and that the interest lists for the LPU and 
the TADD-2 are running consistent with historic 
levels, TAPR will continue to be financially healthy 
through 2009. 

The standard Balance Sheet and P&L reports for 
2008 will be presented at the annual meeting at the 
DCC. 

(Editor's Note: The Treasurer’s Report is 
preliminary and subject to approval by the TAPR 
Board of Directors.) 


Ht 


Write Eary 
AND Write OFTEN 


Packet Status Register (PSR) is looking for a few 
good writers, particularly ham radio operators 
working on the digital side of our hobby, who 
would like to publicize their activities here. 

You don't have to be Vonnegut to contribute to 
PSRand you dont have to use Microsoft Word 
to compose your thoughts. The PSR editorial 
staff can handle just about any text and graphic 
format, so don’t be afraid to submit whatever you 
have to wallou@taprorg. 

‘The deadline for the next issue of PSR is July 
15, so write early and write often. 


HHH 
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A Cheap SDR Loran-C frequency receiver 


Poul-Henning Kamp 


<phk@FreeBSD.org> 
Den Andensidste Viking 
Herluf Trollesvej 3 
DK-4200 Slagelse 
Denmark 


ABSTRACT 


Loran-C radio-navigation signals are broadcast in most of the northern hemisphere 
and offer a solid alternative to GPS disciplined frequency standards, but a regrettable lack 
of affordable receivers means that Loran-C sees very little use for this purpose. 


This article describes a simple, cheap and efficient Loran-C frequency receiver for such 


purposes. 


The three main components of the receiver is a homebuilt loop antenna, the TAPR/N8SUR 
"ClockBlock" and the OliMex.com ADUC-P7026 microcontroller prototype card, for a 


total cost well under $200. 


Figure 1. The complete receiver. 


1. Introduction 


LORAN-C signals, like GPS signals, offer 
wireless access to very stable and calibrated, and 
thus expensive, frequency sources and can there- 
fore be used to steer cheaper oscillators to correct 
frequency. 


Despite stern warnings to the contrary, the world 
today is more or less entirely frequency locked to 
the GPS satellites because that is the most con- 
venient and cheap technology available. 


This little project attempts to show that an equally 
cheap solution can be made based on the Loran-C 
signals, which are broadcast over most of the 
northeren hemishere. 


Compared to GPS, Loran-C signals are just about 
as different as can be, where GPS satellits trans- 
mit their signals with only SOW using 1.2 GHz 
microwaves from 20,000 km above the Earth, 
Loran-C transmits are 200 meter tall steel towers 
which transmit several hundred kW at a frequency 
of 100 kHz. 

This makes Loran-C the perfect backup for GPS, 
since there are almost no overlap between threaths 
to the two systems. 

This paper documents what is clearly a "work in 
progress”, and as time and energy permits, I will 
update both this paper and the source-code that 
goes with it. 

Poul-Henning Kamp 


108 SPRING 200 


2. Hardware 


The entire reason for this article is that a new 
breed of microcontroller combines a real 32bit 
CPU with fast analog/digital converters on a sin- 
gle chip. 

Amongst these chips my fancy took to Analog 
Devices ADuC7026. 


beecdty bea 


HP EE EE 


Figure 2. ADuC7026 block diagram 


For one thing, Analog Devices are experts at the 
analog/digital border, and the specifications on the 
ADC in the ADUC is much tighter than on on 
competing chips. 

Furthermore the ADC offers a very convenient 
fully differential input mode and is not picky 
about the ADC reference voltage. 


There are also some downsides: The chip needs a 
42 MHz clock frequency to run the ADC at 1 
megasample per second, and it only has 8kilo- 
bytes of RAM. 

A very interesting feature is the built in "pro- 
grammable logic array", which makes it possible 
to totally avoid external glue-logic in our case. 
The main market for the ADuC7026 seems to be 
the variable frequency motor control market, and 
certain on-chip facilities are less than well 
thought out and documented. 


One example of this is that the internal timers 
barely have any connection to the outside and that 
trigger conditions have unexplained systematic 
delays. 

As it transpires, there is a way to hook things up 
in a way we can use. 


2.1. The 42MHz clock frequency 


The 42MHz clock frequency is much less of a 
problem than it sounds, thanks to John Acker- 
mann, aka N8UR, we can simply pick up a 
"ClockBlock" from TAPR.org, set the jumpers 
correctly and get 42MHz out of almost any input 
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Hardware 


frequency we care to use. 


Figure 3. TAPR.org ClockBlock card 


It is important to jumper the ClockBlock for 3.3V 

so it does not overdrive the ADuC7026 clock 

input. 

The calculation of correct jumper settings is semi- 

complicated, and by far the easiest way is to use 

the web-page IDT provides: 

http: //timing. idt .com/ 
calculators/ics525inputForm. html 

For your convenience, here are three typical 

jumper settings: 


| F,= | IMHz SMHz 10MHz 
ST) 0 0 0 
SI 0 0 0 
$2 1 1 k 
RO 0 0 0 
RI 0 0 0 
R2 0 0 0 
R3 0 0 0 
R4 O) 0 0 
RS 0 0 0 
R6 ) 0 0 
vo 0 0 1 
vi t 1 0 
v2 0 0 1 
v3 1 0 1 
v4 ) 0 0 
vs 0 1 0 
V6 1 i) 0 
v7 1 0 0 
vs ) 0 0 


Figure 4. ClockBlock jumper settings 


2.2, The ADuC7026 microcontroller 

Being mostly of the software persuation, I find the 
prospect of soldering 80 pin SMD packages 
something to be avoid if reasonably possible. 


Fortunately my favourite Bulgarian electronics 
pusher, OliMex.com, has an ADUC-P7026 proto- 
type card which is perfect for this project. 

In the USA, you can purchase it from Spark- 
fun.com. 


Figure 5. Olimex.com ADUC-P7026 card 


2.3. The Antenna 


I do not really consider the antenna part of the 
project, because I reused a very simple loop 
antenna I built for a previous round of Loran-C 
experiments, based loosely on a diagram I found 
on the website VLF . it, using an AD797 opera- 
tional amplifier. 


Figure 6. $20 homebrew loop antenna 


It goes without saying that I have not spent a lot 

of time on antennas, but roughly speaking, the 

requirements for the antenna input signal are 

something like: 

cA Relatively flat passband from at least 70 to 
130 kHz. 


Not too much signal above 300-400 kHz. 
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O At least 100mV peak-to-peak Loran-C signal. 
O No more than 3V peak-to-peak total output 
signal. 


2.4. The Antenna Balun 

The ADC in the ADUC 7026 microcontroller can 
operate in a number of modes, but by far the most 
convenient for our purposes is the truly differen- 
tial mode. 

A differential CMOS ADC can be conveniently 
driven by a signal transformator of the kind found 
in old ADSL or ISDN equipment. 

One particular nice thing about this approach is 
that a transformator has offset voltage. 


ANTENNA BALUN 


ANT 
x4 ADc2 
RA 27R 
>—— t ' i ahi R3 ss 
* 7 8 C3 
} Al 
SH 2 278 x 
a4 R4+ 
4zen —_ AVDD 
NB: R5_= 5OR x ct, A 


(TRL ratio*2 is I: 
gUS “T2100 AGND 


Figure 7. Antenna Balun Schematic 


To be honest, this balun consists of whatever I 
found in my junkbox and is not optimized in any 
way shape or form. 

The transformer I found is a 1:2 trasnformer from 
an ISDN NT box, where both sides have center 
taps. I feed only one half of the input side so it 
looks like a 1:4 with secondary center-tap. 


Notice that R5 should be chosen as the input 
impedance multiplied by the square of the trans- 
formators ratio: 


RS = Zin > N? 
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Fin [TAPR/NSUR| 42MHz| - ADuC7026 microcontroller 
™| ClockBlock = H 1 } 
' ' 
<| PWM FIQ |ARM7TDMI| | 
1MBz| | cru. |! 
ee Timer! } 

' _ |/-»} PLA[S] | »!Mark signal 
BALUN e<nai ADC Ie 
\ Vref \ PLA[6] |_| 
gbviet J ' 
H \ 

le} DACO | CoStart | PLA[7] }»!Debug signal 

H i 
H i 

(\)-—-J DACI || DAC2 || DAC3 |®\EFC signal 
Soe ' \ 
Oscilloscope Lisscocassassccesesecesce 4 


Figure 8. Blockdiagram 


2.5. Connecting it all 


That is really all for the hardware, now we just 
need to connect it together: 


© Jumper the ClockBlock for 3.3V supply. 


© Connect the ClockBlock output to pin P0.7 on 
the ADUC-P7026. 

Connect your chosen frequency standard to 
the ClockBlock. 

Connect pin P2.7 (PWM1;) to P0.6 (TI) on 
the ADUC-P7026 to establish a connection 
from the PWM output to the Timer] input. 
Connect pin XXX to XXX on the ADUC- 
P7026 to establish a connection from 
DACXXX output to the ADC’s VREF input. 
Connect the differential antenna signal to pin 
ADC1 and ADC2, remember that they both 
must stay between 4,,4 and A,iq at all times. 
Find a suitable power supply for both the 
cards, it should be at least 6VDC to avoid 
drawing power from the USB port, and needs 
to supply about 100mA total. 

Connect your computer to the USB port on 
the ADUC-P7206 card. 


o 


2.6. Available output signals 


The following output signals are available for reg- 
ulation and debug use. 


of the receiver hardware. 


a P2.1 (PLAO[6)) "Mark" output. The raising 
edge marks the start of measurements in the 
A-code GRI interval. The software reports 
how much later than this flank the 3rd zero 
crossing is found. 

P2.2 (PLAO[7]) "Debug" output. This pin 
flips low whenever and as long as the fast 
interrupt routine runs. 

PO.S (ADC ius») This signal is logic high while 
the ADC is performing a conversion 
DACXXX Analog replica of the averaged sig- 
nal. Show this on your oscilloscope, use the 
"mark" output as timebase trigger. 


2.7. Bootloader trigger circuit 


The ADUC-7026 has a very convenient serial 
boot-loader which I use to download the firmware 
to the internal flash memory. 


Unfortunately, there is no way to enter the boot- 
loader from software, it can only be triggered if 
P0.0 is logical low during a reset. 


This little circuit implements a crude mono-table 
timer which can hold P0.0 down for us, while we 
trigger a reset from software, thus saving a trip 
down to the lab. 


APR PSR 


Aduc7@26 
Bootloader peer 
Control 


2.8 
Pat an 
Rt BC547 


GND 


6NO GND GND GNO 
Figure 9. Bootstrap circuit 


3. Loran-C signals 


A Loran-C navigation chain consists of a master 
transmitter and from one to five slave transmit- 
ters, all broadcasting on the same "GRI", Group 
Repetition Interval. 


A Physical transmitting station can participate in 
multiple chains, although for reasons of time-con- 
tention, typically it will only participate in two 
chains. 

In our case, we are interested only in a single 
transmitter in a single chain, typically the nearest 
and strongest signal at the location of reception, 
although there is nothing preventing this construc- 
tion from being used for more ambitious recep- 
tion projects. 


3.1. What the station transmits 


There are many ways to describe the Loran-C sig- 
nal, but for our analysis of how and why this 
receiver works, we will treat it as three compo- 
nents which are convoluted to give the actual on- 
air signal. 


3.1.1. The GRI periodic function 


This is basically a pulse generator which emits a 
pulse every GRI * 10 yseconds, so that for 
instance the Sylt chain with GRI 7499 has a 
period of 0.07499 seconds. 


A very important feature of the GRIs used, is that 
they do not result in periods that are integral mul- 
tiples of 1 msecond, so signals from CW radio 
transmitters will average out, allowing us to use 
essentially no other frequency domain filtering. 


3.1.2. The signal code function 

This is where much of the noise resistance of the 
Loran-C signal structure comes from: in alternat- 
ing GRI periods two different pulse-trains are 
used for each of Master and Slave stations. 
Whenever the GRI periodic function triggers a 
transmission eight or nine pulses will be transmit- 
ted with 1 msec interval, with polarity according 
to the following two figures. 
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Loran-C signals 


Master A-Code 


Master B-Code 
Figure 10. Master Signal Codes 


In addition to the frequency filtering provided by 
the GRI averaging, the non-periodic nature of 
these codes provide significant improvements in 
SIN ratio. 


Slave A-Code 


Slave B-Code 
Figure 11. Slave Signal Codes 


One important feature of both codes is that if the 
A and B interval codes are summed up, the odd 
numbered pulses cancel out leaving only four sig- 
nals, spaced 2 msec apart. 


3.1.3. The Radio Frequency Pulse function 


The final step of the modulation is that each pulse 
is transmitted as 


f@=sin (=) P exp? 


where ¢ has units of useconds. 


| 
i 
ie 28 


Figure 12. Loran pulse 


al i 
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In reality it is only the first part of the pulse-shape 
which is controlled, and in particular, the refer- 
ence point is the 3rd positive zero-crossing of the 
signal 


T 
| 
i 


i 


| 
ae 


Figure 13. Loran pulse front 


When we combine all these components, the 
result is something like this: 


Figure 14. 6731 Lessay chain 


All four stations in the Lessay chain are visible in 
this plot, first the master in Lessay, notice the 
"extra" 9th pulse 2 msec after the 8th, the follows 
the Xray slave in Soustons, the Yankee slave in 
Rugby and finally the Zulu slave on the island of 
Sylt. 

Knowing the exact geographical locations of 
these four transmitters, and the time interval from 
the master transmission until the slaves transmit 
their pulses, we could calculate our geographic 
position from the measured time-intervals 
between these signals. 
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4. Reception 


So, how do we find and lock onto the Loran-C 
signal with only 8 kilobytes of RAM ? 

The brute-force method would be to simply aver- 
age all samples over the FRI interval, that would 
allow us to identify all signals in the chain, tell A 
and B codes apart, locate the 3rd zero-crossing 
and track each signal. 


Unfortunately, this would require up to 
29999 - 10 = 199980 averaging buckets. 


Even if only one bit is used per bucket, this 
amounts to 24.5 kilobytes, more than three times 
the amount of RAM we have available. 


We can get rid of the factor two by averaging only 
over the GRI interval, and provided we rotate 
every other GRI period one millisecond, we can 
still tell the four signal codes apart, at the expense 
of recognizing two energy levels in the pulse pat- 
terns. 


Figure 15. delayed master code average pattern 


Unfortunately, that still requires more than 12 
kilobytes of RAM. 


It is possible to implement both of these versions 
of the algorithm by looking only at a small win- 
dow of the GRI or FRI period at a time. 


Since we realistically need 16, or preferably 32 
bits, per bucket, and at best can use around 6 of 
the 8 kilobytes, this would result in at least 33 
windows, in practice more because of necessary 
overlapping, each of which must be integrated for 
some time, before a signal may be appearant. 

Instead we this narrow scan, single phase 


approach, we opt for a widescan, multi phase 
approach. 
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5. Phase Zero - Locating the strongest signal 


Our first phase must locate the strongest signal in 
the chain, typically the transmitter closest to the 
receiving antenna. 


When averaged over the GRI period, both the 
master and slave codes become four pulses sepa- 
rated by 2ms and our goal is simply to find the 
strongest one of these. 


We do not care at this point, if it is a master or 
slave signal or where the A or B codes are in 
time, the next phase will figure that out. 


Assuming we have 1000 buckets of 32 bits avail- 
able, we can sample the GRI interval every 
100us, and in doing so, we would likely miss 
even very strong signals: We need a sampling fre- 
quency which is not a divisor in the 100 kHz car- 
rier frequency. 

Since our goal is to just find the rough location of 
the strongest signal in the GRI, we are not even 
constrained to sample at rates that divide into the 
Inss pulse spacing, and in fact get better results by 
not doing so. 

Experiments and simulations has shown that if we 
average every 114s, every 114th sample, over 
the GRI period, we can expect to catch at least 
30% of the peak amplitude of the strongest signal, 
and need only 99990 15/114 = 887 buckets. 


5.1. Tek4014 view of phase zero 
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Figure 16. Tek4014 view of phase 0 


This is how the tek4014 window looks in phase 
zero, we see the averaging buckets plotted, and 
the peak signal is marked. 

In the lower left comer of the box is the amplitude 
of the signal in ADC units. 
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Phase zero 


Notice that we do in fact quite clearly see all four 
peaks in the averaged signal, but one of them is 
significantly taller than the other three. 


5.2. Hardware view of phase zero 
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Figure 17. Phase 0 GRI wide signals 
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Trace 1 (black) is the marker output which is used 
to trigger the oscilloscope. 

Trace 2 (green) is the debug output. 

Trace 3 (red) is the DAC] output, this is the same 
data which is plotted in the TEK4014 window. 
Trace 4 (cyan) is the antenna input signal, where 
we can spot the LORAN pulses, thanks to a bit of 
averaging in the oscilloscope. 
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Figure 18. Phase 0 zoomed signals. 


Zooming in, we see a single pulse of the loran 
signal in trace 4. Trace 2 shows the sample dis- 
tance of 114 microseconds and trace 3 shows that 
we almost managed to miss this loran pulse 
entirely between two samples. 
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6. Phase One - Identifying the code 

The peak signal from phase 1 must by definition 
be within 57 samples of the true peak value, of 
one of the four peaks, 2 msec apart, resulting 
from the GRI summing. 

Now we need to find out if this is a master or 
slave station, and to identify which half of the FRI 
has the A and B code. 

If we caught the last of the four peaks, the code 
starts 6 msec ahead of this point, if we caught the 
first of the four peaks, it spans another 9 msec 
after this point, so all in all, there are 16 windows 
we need to check for the precense of pulses. 

We use a repetition frequency of 2 * GRI, some- 
times called "FRI" for Frame Repetition Rate, in 
order to not add the A and B codes together. 
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Figure 19. model code windows 


Using no more than the 876 integers used in 
phase zero, each bucket can get 54 integers and if 
we space them 13 samples apart, each bucket will 
cover 702 useconds, 

Again, 13 is a good number because it will make 
the sample points "wander" over the 100 kHz 
period of the pulses. 
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Figure 20. Tek4014 view of phase 1 
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Phase one 


13 


The tek4014 view shows the 16 windows and 
shows the starting point of the strongest code in 
them. 

Above each window is the energy content of that 
window, in thousands of the maximum found in 
any window. Due to timing effects, it is quite 
likely that one or two windows show numbers 
slightly above 1000. Presently, the energylevels is 
not used for code identification. 

Below the windows are statistics for each of the 
best match for each of the four possible codes. 
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Figure 21. Phase 1 FRI wide signals 


This is the hardware view, showing a full FRI 
interval. 


Notice in trace 4, that we only sample every other 
pulse group. 
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Figure 22. Phase 1 pulse group 


Zooming in on the pulse group, we can see in 
trace 2 how each bucket is sampled and in trace 3 
the resulting waveform. 

Notice that the pulse shape is not fully recovered 
at this point, we only sample every 13 


microseconds, but that is plenty to identify the 
code. 
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Figure 23. Phase 1 single pulse 


Here we are zoomed into just two pulse windows. 
7. Phase Two - Locking on to the signal 


Now we have all the timing informaton we need 
to start integrating the signal for good: we know 
where the maximum is of the first pulse in the 
code, and we know which code it is. 


In this phase we integrate every ADC conversion 
in a 750 (XXX: check code) sample window 
around all the 16 or 18 pulses in the signal, into 
the same single bucket, taking care to invert the 
pulses per the code. 
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Phase two 
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Figure 25. Phase 2 FRI view 
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Figure 26. Phase 2 group view 


Zooming in on the group, we can see th nine 
pulses of the master code being sampled. 

Notice in trace 2, that there is a short burst of 
samples in the space between the 8th and 9th 
pulse, and after the 9th pulse. These samples are 
not technically necessary, but are planned for 
AGC use. 
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Figure 24. Tek4014 view of phase 2 


The tek4014 view shows two views of our inte- 
gration buffer, with the bottom one zoomed in on 
the peak value. 

The various candidates for 3rd zero crossings are 
marked and their statistics given in numeric form 
below. 

The MARK pulses are spaced 2 * GRI apart now, 
and we can see how the signal is sampled twice in 
that period. 
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Figure 27. Phase 2 pulse view 


Finally zooming in on a single pulse, we can see 
the full wave-shape of the Loran-C pulse, includ- 
ing the sky-wave echo. 


The small "pre-echo" is currently unexplained. 
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Figure 28. Finding 3rd zero 


8. Phase Three - Tracking the 3rd zero-cross- 
ing. 

In state 3 we have chosen the zero-crossing we 
want to track, and hold on to it to the best of our 
ability. 

In practice we cannot simply nail a particular 
sample in the GRI window and track that, as the 
signal and the timebase may wander relative to 
each other. 

We always use as tracking sample, the sample that 
is closest to the 3rd zero crossing. 

Nominally the tracking point is sample 75 in the 
window, but we allow it to wander up to 10 sam- 
ples in either direction, before we realign the sam- 
pling window to put it back at sample 75. 

We interpolate the true zero crossing between the 
tracking sample and the two samples that sur- 
round it. 


This is the Tek4014 display in phase three: 
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Figure 29. Tek4014 view of state 3 


The top plot shows the entire sampling window 
and a marker for the tracking sample. 

The bottom plot zooms in and marks both the 
tracking sample and the two neighboring samples. 
Below the plots are the pertinent statistics, includ- 
ing calculated timing of the 3rd zero crossing. 
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Phase three 


The three samples which contribute to the inter- 
polation of the zero crossing, are each filtered 
with a FIR bandpass filter: 


ii 


Figure 30. The FIR filter kernel 
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Figure 31. The FIR filter frequency response 


The exact choice of filter has proven to be 
remarkably unimportant, but it does in fact reduce 
the noise in the three sample points considerably. 


8.1. Interpolating the zero crossing 

Since we want better resolution than the whole 
microseconds the sample rate provides, we need 
to interpolate the true zerocrossing from the sam- 
ple values around it: 
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Figure 32. Interpolating the zero crossing 


If we define Yn as the sample that has the lowest 
absolute value, closest to the 3rd zero-crossing, 
we can determine the range of Tx by solving: 
Loran(x) =—Loran(x + 1s) x € [29us--- 30s] 
Doing so, we find x = 29. 664 us. 

It follows that the Yn sample must be located 


between [29. 664yus--- 30.6645] and conse- 
quently that Tx €[-337ns-- - 664ns]. 
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The analytical solution for finding Tx given Yb, 
Yn and Ya is neither practically realizable with the 
limited amount of CPU and memory we have 
available, nor necessary. 


Instead we calculate the FIR filtered signal at the 
sample nearest the zero-crossing (Fn), and the two 
samples right before (Fb) and right after (Fa) and 
find the ratio: 
Fn 
Fa-Fb 


€ [-0. 2871 --- 0.2679] 


Which we map to the corresponding: 
Tx — [-337ns- --664ns] 


Using a lookup table containing the precomputed 
conversion function: 


Figure 33. Interpolation non-linearity 


The lookup table is built by the script 
pulse_param.tcl and lives in the source file 
pulse_param.c. 

The script automatically reads in the FIR filter 
coefficients from the fir. c source file and deter- 
mines how many entries the lookup table needs to 
have, to provide the specified resolution. 

The table is extended 5% further than necessary 
in both end, to cater for edge effects and noise 
components. 

For a resolution of nanoseconds, the table has 
1111 entries. 


8.2. Serial output data 


age 1324 

avg 10 

track 75 

yb 22376693 
yn 8627207 
ya - 8645207 
Fyb 12611456 
Fyn 2030392 
Fya -10263240 
ps 105 

Phase 0 

frac -185 

off = 185 

# 


Figure 34. Serial data stream 
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Phase three 
9. Performance 


These are some of the first performance data col- 
lected, and they should therefore be taken with all 
sensible precaution. 

The signal is 7499M, 205 km away. 

The exponential averaging factor is 1/1024. 

The timebase is a free-runing FRS10 Rb, the sam- 
ples have been compensated for the randomly 
chosen, but deliberate 7.225e-10 frequency offset. 
A phase measurement is recorded once per sec- 
ond using Timer2, from approx 2008-11-30 11:00 
to 2008-12-01 11:00 UTC. 


Figure 35. Raw offset samples 
The standardeviation of the raw samples is 22.6 
nanoseconds. 
A histogram of the samples show a gaussian dis- 
tribution, which seems well suited to more aggre- 
sive averaging. 
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Figure 36. Offset value interval histogram 
The modified allan variance is plotted below: 


Figure 37. Modified allan Variance 
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10. Software 


The software that goes with this article is writte as 
proof of concept, but is written such that experi- 
mentation and further development is possible. 


10.1. The tricky assembler bits 
The tricky bit of the software is the assembler 
coded fast interrupt routine in the file crt0.S. 


T have written this code so that it does not encode 
any of the high-level logic, but instead mechani- 
cally walks a chain of structures that tell it what to 
measure and when. 

For instance, when we scan the 6731 GRI in 
phase zero, the specification looks like this: 


{ 


mark = 1; 

ptr = ss->ptr; 
polarity = 4"; 
avg = 6; 

ent = 1; 
timer_repeat = 114; 
repeat = 
timer_after = 

next = 


Figure 38. Phase zero specification 


This specification tells the assembler code to 
pulse the MARK pin high when starting this spec- 
ification, then 590 times take one ADC measure- 
ment 114 microseconds apart, then skip 50 
microseconds and start over. 

Since 590 * 114 + 50 = 67310 this will repeatedly 
scan the 6731 GRI signals. 

The polarity member specifies that the samples all 
have positive phase, the avg specifies he exponen- 
tial average time constant and the next member 
tells it what to do next (the same thing). 

With none of few modifications, this scheme 
should also support reception of multiple stations 
in the same chain or even several stations from 
different chains. 


10.2. The high level logic 

The highlevel logic, such as it is, is written in C 
code in the loran@.c file. 

Right now the code is semi-manual, controlled by 
single characters received on the serial/USB port, 
running at 115200 bps. 
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Software 


The following commands are recognized: 


c This will present a menu with Loran-C chains 
from which you can choose one by entering 
the character in the [. 


if 


2] 
a) 
8) 
cl 
D} 
z 
F 
ic} 
cy 
bi 
a] 
K] 
| 
iM 
N) 7980 Southeast U.S. 
B 
S| 
Tr] 
ti) 
iv) 
C) 
x! 
¥ 
2} 
i 


in 
i 
5 


Calcutta 

5930 Canadian Bast Coast 

5980 Russian-American 

5990 Canadian West Coast 
Bombay 


Lessay 
6780 China South Sea 


2 
7030 Saudi Arabia South 
7270 Newfoundland Bast Coast 
7430 China North Sea 


t 
7980 Bistem Russia Chayka 
Gulf of Alaska 


2 
tS 
3 


8930 North West Pacific 
8970 Great Lakes 


Bidi 
9610 South Central U.S. 
9930 Bast Asia 
9940 U.S. West Coast 

9960 Northeast US 

9990 North Pacific 


Please select chain: 
Figure 39. Chain selection menu 


Make TEK4014 plot via serial/USB port 
Zoom in on curves in TEK4014 plots 

Zoom out on curves in TEK4014 plots 
Double the exponential average time constant. 
Half the exponential average time constant. 
Skip to mode 1 

Skip to mode 2 


vecerad 


10.3. TEK4014 plotting 

Before windows and mice, Tektronix produced a 
series of high quality graphical terminals, based 
on the storage-CRT principle. 

These plotting instructions to these terminals 
became the de-facto format for plotting files and 
as a result, the original X11 "xterm" program 
offered, and still does, TEK4014 emulation. 

If you connect to the serial/USB port using the 
X.org xterm program and press °t’, then you will 
be rewarded with a nice plot as seen earlier in this 
paper. 

Tf you use another terminal program which does 
not support TEK4014 plotting commands, you 
will get a blast of random ASCII characters. 


10.4. Utility functions ete. 


The other three C language files contains the code 
to configure the hardware, code to use the 
serial/USB port and a function which calculates 
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basic statistics for an array of integers. 


The files divdi3.c, qdivrem.c and divsi3.S, are 
runtime support files for the GCC compiler. 


10.5. Downloading firmware to the ADuC7026 


You can either use a JTAG gadget or download 
the firmware via the serial/USB port. 


Included in the source-code package is an "aduc" 
program I have written for that purpose. 


10.6. Compiling the firmware 


Cross-compiling is notoriously tricky, but the 
FreeBSD build enviroment makes it quite easy. 

To build a arm cross-compiler + toolchain, simply 
do the following: 


ed /usr/sre 
make toolchain \ 
TARGET=arm TARGET_ARCH=arm 


Figure 40, Building an ARM toolchain on FreeBSD 


The Makefile knows how to find these tools under 
the /usr/obj. 
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Figure 41. Protocol analyser trace 


11. Timer] unexplained 


This is a confession I have no idea exactly how 
Timerl is supposed to work in the ADWC7026 
chip, and the way I have observed it work baffles 
me, but the software manages to work around it. 


The PWM section provides us with a IMHz fre- 
quency which, through an external pin, becomes 
Tim erl ’s clock signal. 

In Fig XXX above this external signal is named 
PWM, which for this experiment, was give avery 
short "on" time. 

The signal named "MARK" is from PLA[5], 
which is configured so the flip-flop is latched by 
Tim erl’s count-complete signal. 


The signal named "ADCBUS" is the ADCbusy 
signal, which indicates that the ADC has started 
performing a conversion, and the ADC is also 
triggered by Timerl ’s count-complete signal. 


So far, so good it looks like Timerl expires 
around the middle of the figure. 


But now look at the "DEBUG" signal, which is 
programatically lowered as the first thing in the 
fast interrupt handler, then pulsed high efter 
Timer! has been reloaded, 


The fast interrupt is also triggered by Timerl’s 
count-complete signal, but this happens well 
before PLA[5] and the ADC receive the signal. 


From the FIQ request is raised until the first 
instruction of the handler is executed takes at least 
five clock cycles, but it can take as much as 13 
(XXX 2 cycles, depending which instruction the 
CPU was executing. 

This is why the "DEBUG" trace shows multiple 
transititions: it isnot deterministic in time. 

But he mystery gets deeper here, because the pre- 
vious PWM signal, presumably the one that 
counds down to one, is a fair bit further shead of 
the DEBUG signal than 5-13 (XXX) clock cycles. 


But it gets more mysterious still: The fast inter- 
rupt handler loads a value into the TILD register 
right before the high pulse on the DEBUG signal 
(tight of the dotted line). 


If the leftmost PWM signal is the 1 count, then 
the central one is the 0 count, it would be reason- 
able to expect that loading N or N-1 into the 
TILD register at this paint, would cause the next 
Timerl count complete event to happen N cycles 


APR PSR #108 SPRING 200 


Tim erl unexplained 


of the PWM signal later, starting with the one in 
the right hand side of the trace. 


Experiementally I have found that I have to load 
N-3 to make that happen. 


For this to make sense, the PWM edge that trig- 
gered the count down to zero must be just outside 
the papers left edge so that the central PWM pulse 
is number 3 in the new counting cycle. 


If that is the case, then the chip applies approxi- 
mately 3us worth of metastability latching on 
Timerls count-complete output, before it reaches 
the ARM7TDMI core and the ADC unit. 


At the 42MHz clock, 3u correspond pretty pre- 
cisely to 128 clock cycles. 


If that is the depth of the latch-line that resolves 
metastability, then I would pressume to think that 
it is sufficient deep. 


Enquiring minds wants to know... 
HHH 


DCC DVDs For Sale 


Did you miss the ARRL and TAPR Digital Communica- 
tions Conference last September? Did you attend, and now 
want a record of some, or all of the talks? 

Good news! ARVN:Amateur Radio//Video News cap- 
tured every minute of the presentations on video, and now 
offers them on six DVDs in high quality video and audio. 

The presentations are grouped on the six DVDs in as 
logical an order as possible. You can buy the DVDs indi- 
vidually ($15/each, plus $3 shipping, $4 international). 
Or get whole set for $75 (plus $3/$4 shipping). The com- 
plete set saves you $30 over buying them individually, after 
shipping cost. You can order by credit card or PayPal at 
www.arvideonews.com. 

HHH 
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